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ABSTRACT 

We investigate the conditions required for the production of electron-positron 
pairs above a pulsar polar cap (PC) and the influence of pair production on the 
energetics of the primary particle acceleration. Assuming space-charge limited flow 
acceleration including the inertial frame-dragging effect, we allow both one-photon and 
two-photon pair production by either curvature radiation (CR) photons or photons 
resulting from inverse-Compton scattering of thermal photons from the PC by primary 
electrons. We find that, while only the younger pulsars can produce pairs through CR, 
nearly all known radio pulsars are capable of producing pairs through non-resonant 
inverse-Compton scatterings. The effect of the neutron star equations of state on 
the pair death lines is explored. We show that pair production is facilitated in more 
compact stars and more massive stars. Therefore accretion of mass by pulsars in binary 
systems may allow pair production in most of the millisecond pulsar population. We 
also find that two-photon pair production may be important in millisecond pulsars 
if their surface temperatures are above ~ three million degrees K. Pulsars that 
produce pairs through CR will have their primary acceleration limited by the effect of 
screening of the electric field. In this regime, the high-energy luminosity should follow 
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a Lhe oc E rot dependence. The acceleration voltage drop in pulsars that produce pairs 
only through inverse-Compton emission will not be limited by electric field screening. 
In this regime, the high-energy luminosity should follow a Lhe oc E ro t dependence. 
Thus, older pulsars will have significantly lower 7-ray luminosity. 

Subject headings: pulsars: general — radiation mechanisms: nonthermal — relativity 
— stars: neutron — 7 -rays: stars 
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1. INTRODUCTION 

The acceleration of particles and the production of electron-positron pairs are widely 
considered to be two critical elements necessary for generating radiation in rotation-powered 
pulsars. In polar cap models (e.g. Arons & Sharlemann 1979, Daugherty & Harding 1996), 
acceleration occurs in the region of open field near the magnetic poles and 7-rays from curvature 
and inverse Compton radiation produce pairs primarily by one-photon pair production in the 
strong magnetic field. These pairs may screen the accelerating electric field through the trapping 
and reversal of one sign of charge, and may be required for the coherent radio emission process. 
In outer gap accelerators (e.g. Cheng, Ho & Ruderman 1986), a vacuum gap develops along the 
null charge surface and pairs are required to provide current flow through the gap which can then 
operate as a stable accelerator. 

In this paper we discuss plausible regimes of pair formation above the pulsar polar cap 
(PC), including the energetics of relativistic particles and 7-rays that cause and accompany these 
regimes. We treat the acceleration of particles within the framework of an approach elaborated by 
Harding & Muslimov (1998, hereafter HM98), Harding & Muslimov (2001, hereafter HM01), and 
Harding & Muslimov (2002, hereafter HM02) which combines rough analytic estimates and simple 
practical formulae with detailed numerical calculations. As an underlying PC acceleration model, 
we employ the general-relativistic version of a space-charge limited flow model developed earlier 
by Muslimov & Tsygan (1992, hereafter MT92) and advanced in a number of important aspects 
by HM98, HM01, and HM02. The main focus of our present study is the physical condition 
for pair formation and how this condition translates into a theoretical pair death line for the 
observed radio pulsars. This paper is a logical epilogue of our previous studies (see HM01 and 
HM02) where we calculated the parameters of pair-formation fronts (PFF), including the flux of 
returning positrons, calculated X-ray luminosities due to PC heating, estimated luminosity of 
the primary beam, and revised derivation of pulsar death lines. In our calculations of PFFs we 
employed the standard mechanism of magnetic pair-production by high-energy photons, generated 
via curvature radiation (CR) and/or inverse Compton scattering (ICS). In HM02 we calculated, 
both analytically and numerically, the theoretical pair death lines based on the abovementioned 
regimes of pair formation. However, in HM02 we presented the results of our calculation of pulsar 
death lines only for a canonical neutron star (NS) with the mass 1.4 Mq and radius 10 km, even 
though we pointed out that the effect of deviation of NS mass and radius from their canonical 
values might be important for our calculations. 

In the present study we extend our previous analysis (HM02) of pulsar death lines to explicitly 
incorporate the effect of different NS mass and radius and, for the short-period (millisecond) 
pulsars, to include the possibility of two-photon pair formation. We must emphasize that the 
effect of bigger NS mass is especially worth considering for the millisecond (ms) pulsars, which are 
believed to be descendants of accreting NSs in low-mass binary systems. It is remarkable, that 
our present study suggests that the ms pulsars do favour the bigger NS masses which seems to be 
consistent with their standard evolutionary scenario. This effect is associated with the dominance 
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of the relativistic frame-dragging term in the accelerating voltage drop which is a unique feature 
of the electrodynamic model of MT92. The frame-dragging component of the electric potential 
(field) is proportional to the general-relativistic parameter k = el /MR 2 , where e = r g /R, r g is the 
gravitational radius of a NS of mass M, and / and R are the stellar moment of inertia and mass, 
respectively. Another important aspect of our previous and present studies is the derivation of a 
theoretical relationship between the pulsar's 7-ray luminosity and its spin-down power/luminosity. 
In this paper we discuss such a relationship in the context of the available and forthcoming 7-ray 
observations of pulsars. 

The paper is organized as follows. In §2 we discuss the determination of pair death lines in 
pulsars. First, we outline the basic definition and main assumptions behind the death-line concept 
(§2.1). Second, we discuss the revised analytic approach in the derivation of theoretical death lines 
(§2.2), and then, in §2.3, we discuss our numerical calculation of death lines. In §2.3.1 we calculate 
the CR and ICS death lines for the NS models with a canonical mass 1.4 Mq and with three 
different equations of state to illustrate the effect of compactness on our death line calculations. 
In §§2.3.2, 2.3.3 we focus on the death line calculations for the short-period (millisecond) pulsars: 
we discuss how the change in mass and radius of the underlying NS model affects the pulsar 
death lines, and present the results of our numerical death line calculation (§2.3.2); we incorporate 
the effect of two-photon pair production for the PC temperatures 1-5 xlO 6 K and present the 
corresponding numerical death lines (§2.3.3). In §2.3.2 and 2.3.3 we illustrate separately how 
the mass of a NS and two-photon pair production, respectively, may affect the death lines for 
ms pulsars. In §3 we discuss energetics of the acceleration of primary electrons and present the 
theoretical relationship between the 7-ray luminosity and spin-down power of a pulsar. Finally, in 
§4 we summarize our main results and discuss their most important implications for pulsars. 

2. Death line determination in PSRs 

2.1. General overview and definitions 

Since the very early attempts to relate the apparent absence of radio pulsars with long periods 
in the P—P diagram with the manifestation of the effect of electron-positron pair formation as a 
condition for their operation, it proved instructive to introduce the term 'death line' to separate 
the domain favouring pair formation from the domain where it would be prohibited (see Sturrock 
1971; Ruderman & Sutherland 1975; Chen & Ruderman 1993, hereafter CR93; and references 
quoted therein). Soon, it became almost a common practice for any theoretical study on radio 
pulsars to produce the resulting death lines. Furthermore, some theories developed the idea that 
on the P-P map it is a death valley (see e.g. CR93) rather than a death line that separates radio 
active from radio quiet pulsars. It is important that during the past decade the number of new 
radio pulsars with a wide range of parameters dramatically increased, which boosted the various 
pulsar population studies. In light of the recent extensive radio pulsar surveys (e.g. Manchester et 
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al. 2001), pulsar population studies, and multi- frequency (from radio to 7-ray) pulsar observations, 
it seems timely to get back to the basic concept of a pulsar death line. 

The standard definition of a pulsar death line implies that pulsar radio emission turns off if 
the energetics of accelerated particles drops below the minimum required for electron-positron pair 
production. Thus, the standard definition of a death line implicitly assumes that pair-formation 
is a necessary condition for pulsar radio emission, and that pulsars become radio quiet after 
crossing the death lines during their evolution from left to right in the P-P diagram. Obviously, 
this basic condition may not be sufficient (see e.g. Hibschman & Arons 2001, hereafter HA01; 
for a most recent study where the sufficient condition assumed was that of pair production with 
high enough multiplicity to compeletely screen the parallel electric field), and functioning of radio 
pulsars may imply far more complex physical conditions (see Usov 2002 for a recent review). 
However, numerous theoretical attempts to produce satisfactory death lines implying even the 
basic necessary condition meet certain challenges. For example, a number of observed ms radio 
pulsars fall below most theoretical death lines. Also, many normal radio pulsars tend to be below 
their death lines based on CR pair-formation. For this reason, and to minimize the underlying 
model assumptions, our previous (see HM02) and present analyses of pulsar death lines are based 
on the minimal requirement regarding pair formation. Note that in all our studies we assume a 
centered-dipole magnetic field of a NS. 

In this paper we illustrate how the spread of NS masses and radii may affect the theoretical 
death lines for ordinary and, most importantly, for ms radio pulsars. Needless to say, compactness 
of a NS is an important parameter in our calculations of particle acceleration (mostly because the 
accelerating electric field is of essentially general-relativistic origin) and pair formation (because of 
a bigger deflection of photon trajectories in the gravitational field of a more compact NS). So, the 
detailed analysis of the effects of stellar compactness on the results of such calculations would be 
quite interesting by itself. However, this effect is worthy of special consideration in the case of ms 
pulsars. The main reason is that the latter are believed to descend from accreting NSs in low-mass 
binary systems and may represent post-accreting relatively massive NSs. As will be demonstrated 
in Section 2.3.2, the increase in NS mass considerably facilitates pair formation in short-period 
pulsars, thus pushing the corresponding death lines down to or below the observed (P,P) values 
for ms pulsars. 

We also consider the effect of two-photon pair production, where 7-ray ICS photons interact 
with thermal X-ray photons from the NS surface to produce an electron-positron pair. This 
process is of primary importance in outer-gap models (e.g. Romani 1996; Zhang & Cheng 1997) 
for pulsar high-energy emission. Zhang & Qiao (1998) investigated the importance of this process 
above PCs of normal pulsars, noting that unrealistically high PC temperatures were required. As 
we will show in this paper, two-photon pairs can be produced more easily above PCs of ms pulsars 
because the PC size is larger. We will compute the death line for two-photon pairs as a function 
of PC temperature. 
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2.2. Analytic death lines 

In this Section we generalize the analytic expressions for the death lines derived by Harding 
& Muslimov (2002; hereafter HM02) to explicitly include the effect of different mass, radius, and 
moment of inertia of a NS. 

To formulate the analytic death line condition we need to know the distribution of voltage 
drop in the pulsar's PC acceleration region. For a given distribution of voltage we can calculate 
the characteristic Lorentz factor of a primary electron accelerating above the PC as a function of 
altitude z and pulsar parameters B and P, j acc (z, B, P). The accelerated electron generates (CR 
and/or ICS) photons that may pair produce if the condition for the corresponding pair- formation 
process (in most cases magnetic pair production) is satisfied. In our previous papers (see e.g. 
HM98, HM01, and HM02) we have demonstrated that use of the pair-formation condition allows 
us to determine the pair-formation altitude as a function of pulsar parameters B and P alone. 
Then, the Lorentz factor of a primary electron evaluated at the pair-formation altitude determines 
a minimum Lorentz- fact or an electron should achieve to generate a pair-producing photon, 
7min(-B, P)- Thus, the death line condition would require that the Lorentz-factor of an accelerating 
primary electron is equal to 7 m i n . In our numerical calculation of death lines we can easily keep 
track of the fulfillment of this requirement and plot the corresponding points in the P—P diagram 
that constitute the death line (or rather death curve). However, an analytic derivation of the 
death line needs an additional independent relationship between the pair-formation altitude z and 
pulsar parameters B (or P) and P. In our previous paper (HM02) we have demonstrated that 7acc 
can be expressed as a function of B , P and an additional parameter, the efficiency of converting 
pulsar spin-down power into the luminosity of the primary beam, 

/prim = ^prim/ Eroti 

where L pr \ m is the luminosity of the primary electron beam, and E rot is the pulsar spin-down 
power (= fi 4 I? 2 .R 6 /6c 3 /(l) 2 , where Bq/ f{\) is the surface value of the magnetic field strength 
corrected for the general- relativistic red shift; all other quantities have their usual meaning and 
will be defined below; see also HM02 for details). 

For the typical radio pulsar parameters P and P and for most obliquities, excluding the 
pure orthogonal case, the dominant term in the expressions for the electrostatic potential and 
electric field in the general-relativistic version of the space-charge limited flow model (MT92) is 
proportional to parameter k. In this paper we shall use the following general expression for the 
parameter k to include its explicit dependence on NS radius and moment of inertia 

k = el /MR 2 = 0.15 x hs/Rl (2) 

where e is a NS compactness parameter, Z45 = J/10 45 g-cm 2 , R§ = R/10 6 cm; M, R and / are NS 
mass, radius and moment of inertia, respectively. 
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Then, the explicit expression for 7 acc can be written as (cf HM02, equation [51]) 

7acc = W 7 f pTim B 12 p- 2 , (3) 

where 

/prim — /prim ^0.15 

and the z dependence of 7 acc in implicit in / pr i m . Here / pr i m is the efficiency of converting 
spin-down luminosity into the luminosity of the primary beam, calculated for a canonical NS mass 
M = 1.4 Mq, radius R = 10 6 cm, and moment of inertia I = 10 g-cm 2 ; and Ato.is = k/0.15; 
B\2 = -Bo/10 12 G, -Bo is the surface value of the magnetic field strength; P is the pulsar spin period 
in seconds. Since the efficiency / pr im depends on the NS radius and moment of inertia, we simply 
normalize it by f pr \ m to make more illustrative the comparison with our calculations performed for 
a canonical NS model. Note that in formula (||) the parameter / pr i m is assumed to be independent 
of B and P, but scales with k and Rq in the same way the 7 acc (or the corresponding acceleration 
potential drop, see eq. [49] in Paper I) does. 

A major advantage of the above formula for j acc is that it does not discriminate between 
unsaturated and saturated regimes of acceleration of primary electrons (see HM02). It implies 
that the Lorentz factor of an accelerating electron is merely proportional to the maximum voltage 
drop above PC, with the coefficient of proportionality being the bulk efficiency of the pulsar 
accelerator, / prrm . 

Now we can write the pair-formation condition as 

Tacc ^ 7miri) (5) 

where 7 acc is given by formula (||). 

The expressions for 7 m i n for different underlying mechanisms for pair-producing photons are 
the same as derived by HM02 but including the explicit dependence on R and k, and read 

Curvature radiation 



{ 3 4 R 5/1 V /7 PV14R-V7 p<p {CR) 

7min I L2 ^1/4 p _ 1/4 p>ptP") t W 



Resonant ICS 



7min 1 0.2 RI'\\%P-^B^' 2 P Z P!; R) , 1 j 



Non-resonant ICS 



-7- 



f R 1 / 2 ^ 1 / 3 P-1/6R-V3 p <p (NR) 

7min 1 0.6 Pf P^P^), W 

where 

pj Cfi) = 0.1 Bfi 9 , (9) 

pj fl) = 0.1 B% 7 , (10) 

P*^ = 0.4 < 7 (11) 

are the critical spin periods (see eqs [4]- [6] in HM02) in the criterion defining the unsaturated 
(upper row) and saturated (lower row) regimes of primary electron acceleration. 

Now, let us use formula (||) in criterion (|5|) to get explicit conditions representing the death 
lines. Note that, following the reasoning of HM02, in criterion (|5|) we should evaluate 7 acc at 
/prim = /prim; w here /p^; m is the minimum pulsar efficiency needed for pair formation. It is this 
minimum or threshold value of / pr i m that determines the pulsar death line condition. We must 
note that at the pulsar death line, in the case of ICS, the /p^; m defines the voltage drop at the 
PFF not the final energy of the primary beam (see Section 3). For each of the mechanisms of 
generation of pair-producing photons, we discuss in this paper, the resultant analytic death line 
(or rather parameter space with allowed pair formation) in the P—P diagram reads 

Curvature radiation 

U P > I f Ig P ~ I !g fA ~ A f R) ( R > J ) " 146 P ~ P * CR) > n 2 ) 
" 1 | lgP-21g/- i m -A^ ) (P,/)-15.4 PZ P ( CR \ 

where a[ CR \r,I) = 1.5(lg/ 45 - 6.31gP 6 ), and A^ R) (R,I) = 1.5(lg/ 45 - 7.31gP 6 ). 
Resonant ICS 

V p> / I lgP-lg^-|A{ JC5) (i?,J)-16.6 P<rt R \ 
Non-resonant ICS 

lg P > / i 1^ ^ — 1 1^ fprhn ~ a[ ICS \r, I) — 18.6 P<P^ R \ 
~ \ | lg P - 2 lg f™Z - A|; C5) (P, /) - 20.0 P ^ Pi^ , 

where a[ ICS ^ = (lg/45 — 7.51gPe) an d A^ S ^ = (lg/45 — 10.5 lgi?6)- in the above death-line 
conditions we used formula (|2|) for k. 

The expressions (0)-(14) differ from the similar expressions presented in HM02 (see eqs 
[52]-[54]) by the terms A's which take into account the deviation of NS radius and moment of 
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inertia from the canonical values of 10 6 cm and 10 45 g-cm 2 , respectively. Thus, for canonical NS 
parameters, the above expressions translate into expressions [52]-[54] of HM02. One can see from 
expressions for A's that the more compact the NS is, the lower the death line moves. 

The analytic expressions above for the ICS pair death lines differ significantly from those 
derived by Zhang et al. (2000). The reasons for these differences were discussed in detail in HM02. 



2.3. Numerical death lines 



2.3.1. Effect of NS equation of state 



The details of the method we use to numerically compute pair death lines can be found 
in HM02. Briefly, we keep track of the total distance, the sum of the acceleration length and 
the pair production attenuation length of either CR or ICS radiated photons, as the primary 
electron is accelerating. The minimum of this total distance is assumed to determine the height 
of the pair formation front. As the value of surface magnetic field decreases for a given pulsar 
period, the PFF moves to higher altitude. Performing this calculation for a range of pulsar 
periods, we find the value of surface magnetic field below which a PFF cannot form within the 
pulsar magnetosphere. This occurs because both the required acceleration length and the pair 
attenuation length become very large. The result is a line in P-Bq space which we identify as the 
death line for pair production by photons of a given radiation type. In order to compare death 
lines for different equations of state (EOS) with the observed pulsar population we must transform 
the calculated lines to P-P space using the magneto-dipole spin-down relation 

1/2 



Bn 



3c 3 IPP 

2vr 2 i? 6 



giving 



P = 2.43 x 10" ib -A 




s • s 



(15) 



(16) 



In our numerical death line calculations we employ three most representative NS models 
standardly used in the calculations of thermal evolution of NSs (see e.g. Table 1 in Umeda et. 
al 1993 and references quoted therein) plus a strange star model (see e.g. Glendenning, 1997). 
The NS models correspond to a star with the baryon mass 1.4 Mq, and different radii and 
moments of inertia: Rq = 1.6 and Z45 = 2.2 (Pandharipande-Pines-Smith'76 model); R§ = 1.1 
and J45 = 1.2 (Friedman-Pandharipande'81 model without pion condensate); and Rq = 0.8 and 
I45 = 0.6 (Baym-Pethick-Sutherland'71 model). The strange star model has a mass 1.4 Mq, 
radius Rq = 0.7, and moment of inertia I45 = 0.7. In Figure 1 we show the death line calculations 
based on these models. Note that the stellar models were produced for a non-rotating star. Thus, 
strictly speaking, our calculations of death lines shown in Figure 1 are not very accurate for the 
ms pulsars. However, as it will be discussed in the next Section our death line calculations based 
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on a non-rotating NS (and perhaps strange star) model may still be satisfactory even in the ms 
range. We should also mention that the only purpose of our inclusion of a strange star model is to 
demonstrate the effect of extreme stellar compactness on our death line calculations. The surface 
physics of a strange star may be significantly different from that of a NS, and in this paper we 
refrain from any speculation on this issue. It was suggested although that some radio pulsars could 
well be strange stars rather than NSs (see e.g. Xu, Qiao & Zhang 1999; Kapoor & Shukre 2001). 

As was suggested by our analytic expressions in Section 2.2, EOS with smaller radii will 
move the death line lower, allowing a greater number of pulsars to produce pairs. Pair production 
is thus facilitated in more compact stars with bigger re's and having softer EOS or even having 
strange matter EOS. We demonstrated this effect by employing the BPS NS model and more or 
less typical strange star model available in the literature. The death line corresponding to the 
strange star model (having largest compactness parameter) is the lowest one of those shown in 
Figure 1. Note, that the ICS pair death lines are more strongly affected by a change in NS radius 
than are the CR pair death lines. In this paper we have computed ICS death lines for only one 
PC temperature of 10 6 K in order to compare the effect of EOS. In HM02, we showed that PC 
temperature has only a small effect on ICS pair death lines for normal pulsars and is much less 
significant than the effect of different NS EOSs. For a PC temperature of 5 x 10 6 K and canonical 
NS model, the ICS pair death line lies slightly below the BPS model death line shown in Figure 1. 

2.3.2. Effect of NS mass for death lines in ms PSRs 

In our calculations of death lines for pulsars with the periods in the range of 0.001-0.1 s we use 
rapidly rotating NS models produced by Friedman, Ipser & Parker (1986). In Figure 2 we present 
our calculated death lines for NSs with the gravitational masses 1.26, 1.97, and 2.64 Mq. This 
particular sequence of rotating NS models is calculated by employing the Pandharipande-Smith'75 
EOS and corresponds to the NS spin period of ~ 2 ms. The NS radii and moments of inertia for 
this sequence are, respectively, Rq = 1.59 and I45 = 2.28, Rq = 1.59 and J45 = 3.9, and Rq = 1.47 
and I45 = 5.18. Note that we used the same sequence of models (implying the NS spin period « 2 
ms) to calculate the death lines for the whole range of spin periods up to 0.1 s. In fact, for the 
spin periods ~ 10 ms the rotating NS models practically converge with the non-rotating models. 
Note also that, for these three models the relative differences between non-rotating and rotating 
sequences in terms of, respectively, the gravitational mass, radius, and moment of inertia, are 3%, 
10%, and 8% (for 1.26 M Q model); 1.5%, 5%, and 5% (for 1.97 M Q model); and 6%, 8%, and 4% 
(for 2.64 Mq model). The differences between rotating and non-rotating models of this magnitude 
are more or less typical for other sequences of models based on a reasonable EOS. Thus, the effect 
of rotation, by itself, is not very important for the death line calculations, and we can justifiably 
use this particular sequence of models for our death line calculations for the period range under 
consideration. What may actually be important for the death line calculations in ms pulsars is the 
mass of a NS. We find that our numerical calculations of death lines in ms pulsars favour more 
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massive NS models, which is consistent with the hypothesis that the ms pulsars descend from 
accreting NSs in low-mass binaries. Figure 2 illustrates the effect of NS mass on the death line 
calculations for the ms pulsars. It shows that the increase in NS (gravitational) mass by 0.6-0.7 
Mq moves the death line down by a factor of a few or more (see also formula [jl4| , saturated case) . 
Even though the mass of 2.64 Mq we use in our calculations may seem to be rather large, the 
main result of our calculations shown in Figure 2 is that the increase in the NS baryon mass by 
~ 30 — 60% may significantly facilitate the process of pair formation above the PC in a ms pulsar. 
This effect may account for the fact that many ms pulsars tend to scatter below the theoretical 
death lines. The fact that the ms pulsars might be more massive NSs processed in binary systems 
could naturally explain this effect. 



2.3.3. Effect of two-photon pair formation 

We also investigate the effect of two-photon pair production on the pair death line. The 
process we consider is that of ICS photons interacting with soft X-ray photons from the hot PC, 
drawing from the same pool of thermal photons that are responsible for creating the ICS photon 
spectrum. In PC models, two-photon pair production has traditionally not been considered 
important in comparison to one-photon pair production. Zhang & Qiao (1998) noted that 
two-photon pair production could be important in PC models if the temperature of the PC was 
high enough (> 4 x 10 6 K for P = 0.1 s). Zhang (2001) estimated the photon attenuation length 
for two-photon pair production above a hot PC of radius Rt to be 

£ 27 ~4.7x 10 5 T 6 - 3 [ 5 (z)S(e)]' 1 ,cm (17) 

where 

g(z) = 0.27-0.507/x c + 0.237/x 2 , (18) 

z 

He = j 

\l* 2 + 4 

where z t = Rt/R, T 6 = T/10 6 K, £(e) = (vr 2 /3) ln(0.1179e)/(9e), 9 = T/mc 2 , and e is the photon 
energy in units of mc 2 . Near the NS surface [i c ~ 0, and at threshold, 

2 . . 

eth = 7^ r, (19) 

9(1 - Mc) 

where e ~ 1/9, S(e) ~ 1 so that £27 — 1-7 x 10 6 T 6 _3 cm. For surface temperatures T§ ~ 1, 
the photon attenuation length, ^27; is much larger than the acceleration length required for 
the electron to radiate an ICS photon above threshold. Therefore, £27 sets the distance to the 
two-photon PFF. Since the soft photon density declines with height above the surface on a scale 
roughly equal to Rt, a reasonable criterion for two-photon pair production is then £27 < Rt- Since 
Rt = r pc = (£IR/ 'c) 1 / 2 R for a heated PC, this condition becomes 

T 6 ~ 1.6(P/1 ms) 1/6 Rl /2 (20) 
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For normal pulsars, Tq > 4 — 5 is required, which is unrealistically high, but for ms pulsars the 
temperature required for significant two-photon pair production is in the range detected for some 
ms pulsars. It is clear that the advantage ms pulsars hold over normal pulsars in the facilitation 
of two-photon pair production is a large PC size, which allows both larger angles between the ICS 
7-rays and the thermal PC photons thus lowering the threshold energy for producing a pair, and 
an increase in the scale length over which the soft photon density decays. The primary electrons 
therefore can reach the energies needed to radiate photons at threshold in a shorter distance. 

In order to compute the two-photon pair death lines numerically, we need an expression for 
the rate of pair production of a high energy photon of energy e: 

i?2 7 (e, 0) = cj d<j) J dfHi J de s a 2 -y(w) n s (e s , fi s ) (1 - /i ti ), (21) 

where [in is the cosine of the polar angle between the propagation direction of the two photons, 
and a2j(w) is the cross section, 

2 

a 2 -y(w) = \(2w A + 2w 2 - 1) ln(u> + Vw 2 - 1) - w(w 2 + l)y/w 2 - lj (22) 

in the center of momentum frame in terms of the variable, 

w= [ee s (l- fxu)^} 1 / 2 . (23) 

The above cross section does not take into account the effect of the strong magnetic field near the 
NS surface. Although these effects may be significant in the highest pulsar fields, the magnetic 
two-photon pair production cross section is very complicated (Kozlenkov & Mitrofanov 1986) and, 
since the process will only be important for ms pulsars having low fields, we will not consider these 
effects here. The above (field-free) cross section may be simplified in two limits, near threshold 
and for large w (Svensson 1982): 



J irr 2 y/w 2 - 1 w ~ 1 , , 

<T2yW ~ | (7rr2>6) [2H 2w) - 1] W » 1, 1 ' 

where r e is the classical electron radius. We choose the coordinate system so that the z-axis is 
along the magnetic pole. To simplify the geometry of the calculation, we assume that the 7-ray 
travels along the positive z-axis, and assume that the soft photons are uniformly radiated over the 
PC. There is thus azimuthal symmetry about the magnetic pole and the polar angle fin ranges 



from to fj, c , where [i c is given in equation (18). The thermal photons from the PC are described 
by the blackbody distribution, 

n s (6,) = (l-Mc)| [exp(ea £ ; e) _ ir (25) 

where A c is the electron Compton wavelength. Changing variables from im to w, and using the 
expressions for (T2^(w) defined by equation (p4|), the expression for the rate in equation (^l|) 
becomes 

R 2l {e) ~ IQ^rl ~ 2 f°° de s S( ^) ( 26 ) 
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where 



[|K 2 -i) 5/2 



l(w 2 s - l) 3 / 2 ] w s ~ 1 



T,(w s ) 



(27) 



2.39 - 



2 In 2w e 



1 



W s > 1 



and 



u; s = max[l, e e s 



(l-A*c)/2]. 



(28) 



Equation (|2(]) is then integrated numerically to obtain the two-photon pair production and 



attenuation length. 

As in the case of the one-photon PFF calculation, we minimize the sum of the acceleration 
length and the pair production attenuation length of ICS radiated photons, as the primary electron 
is accelerating. Performing this calculation for a range of pulsar periods, we find the value of 
surface magnetic field below which a PFF cannot form within the pulsar magnetosphere. Figure 
3 shows the computed pair death lines in P-P space that include the possibility of two-photon 
pair production for different values of the PC surface temperature. We display three cases for 
illustration: 1) death lines for one-photon pairs only 2) death lines for two-photon pairs only 
and 3) death lines for one-photon and two-photon pairs. All cases assume a canonical NS model 
with M = 1.4Mq, J45 = 1 and R = 10 km. It is apparent that two-photon pair production is 
not important at all for any of the known radio pulsar population unless the PC temperature 
Tq ~ 3. The position of the two-photon death line is sensitively dependent on PC temperature 
for 3.0<T6<5.0 and then saturates at about Tq ~ 5.0, reflecting the effect of the two-photon pair 
threshold. For Tq < 3.0, the ICS photons never reach pair threshold during the particle acceleration. 
For 3.0 <Tq< 5.0, the photons are pair producing near threshold where the cross section is sharply 
rising, and for Tq ~ 5.0 the photons are pair producing above threshold where the cross section is 
decreasing. The two-photon death lines curve upward to become almost vertical with increasing 
P because for longer periods, particles must accelerate to high altitudes to reach pair threshold 
where the thermal photon density is declining. Thus, as we had noted previously, two-photon 
pair production is only important in short-period and ms pulsars. Since young, short-period 
pulsars with high magnetic fields do not have detected PC surface temperatures as high as Tq ~ 3, 
two-photon pairs are effectively not important for any but ms pulsars. The combined one-photon 
plus two-photon death lines blend into the one-photon death lines as one-photon pairs dominate 
at higher fields and longer periods. 

HM02 found that substantial PC heating by trapped positrons returning from an ICS pair 
front can occur in ms pulsars if PC temperatures exceed Tq ~ 1. In order for ms pulsars to sustain 
these high temperatures through PC heating, the heated area must be much smaller than the area 
of the standard PC, which is A pc = TrR 2 (ClR/c). This is in fact consistent with the non-uniform 
heating distribution found by HM02. However, for the PC temperatures Tq > 3 needed for 
two-photon pair production the question of the stability of two-photon PFFs must be addressed. 
Positrons returning from the PFF will radiate ICS photons which can produce two-photon pairs 
in a relatively small distance because the pair production threshold (see eq. |l9|]) is much lower 
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for head-on collisions. Creation of enough two-photon pairs by the returning positrons at high 
altitudes could disrupt the acceleration of the primary electrons. Investigation of this effect will 
require inclusion of full angular dependence of the pair production rate and will be considered in 
a future paper. 



3. Acceleration and 7— ray luminosity 

Establishing the regimes of pair formation above pulsar PCs is not only important for 
understanding the behavior of the radio emission, but also allows us to predict regimes of particle 
acceleration and thus high-energy emission since the acceleration of the primary particles may 
be limited by screening at a PFF. HM02 found that CR pairs are very effective in screening the 
electric field at the PFF, whereas ICS pairs are less effective and may only screen the electric field 
above the PFF in some cases. They found that when ICS screening does occur, it only screens 
the field locally but will not screen at higher altitudes. Thus ICS pairs may retard but do not 
ultimately limit acceleration of the primary electrons, which may then also produce CR pairs 
at higher altitude. In fact, the luminosity of the thermal X-rays from a hot PC detected from 
PSR B1929+10, a 3 Myr old pulsar where a detectable cooling component is not expected, would 
require and is consistent with heating by positrons produced at a CR pair front (HM01) since the 
heating by positrons produced only at a ICS pair front would not be detectable (HM02). 

The luminosity of the primary electron beam in the PC pulsar model can be calculated as 

iprim = CiC / \p e \$dS, (29) 



where 

w = wl ( (30) 

is the value of an electron charge density calculated at cos x ~ 1 (where % ls the pulsar obliquity) , 
5>(£,£, <f>) is the electric potential, and 

Q E>3 

dS = ——rfidid<P (31) 
is the element of a spherical surface cut by the last open field lines at the radial distance r (= Rrj). 



Here a = J 1 — r g /R , r g is the gravitational radius of a NS; c is the velocity of light; z is the 
altitude above the PC in units of stellar radius; £ is the magnetic colatitude of a field line scaled 
by the magnetic colatitude of the last open field line; and 4> is the magnetic azimuthal angle. 



In our previous papers (see e.g. HM02) we calculated L pr i m using in formula ( P9| ) the 
expression for the electric potential evaluated at the relatively smaller altitudes (both for the 
unsaturated and saturated regimes of acceleration of primaries) where the bulk of the CR 
pair-formation and electric field screening occur. In the regime where CR pairs are created, i.e. 
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above the CR death line, the luminosity of the primary beam is therefore set by the CR pair front. 
HM02 derived the following expressions for the luminosity of the primary electron beam based on 
the altitude of the CR pair front, 

L (CR) _ 10 i6 f __ /s W 2 £1/2 { P 1/U B^ 7 P<Pi CR \ 

L prim ~ iU ( er g/ S J ^rot i q 3 p-1/4 p > p(CR) ^ 

In the case where there are no pairs produced by CR (and therefore no electric field screening) the 
most appropriate expression for the electric potential in this case is (see eq. [24] in HM01, and eq. 
[13] in HM98) 

flR , 9N f \z r pc /R < z < 1, 

osx f(l)c ^ ' \ l(l-^r) »7-l> 1, 

where r/ = r/i?. This formula applies for the altitudes much greater than the PC radius and 
corresponds to the saturated and unscreened regime of acceleration of primaries. However, in some 
cases the acceleration of primary electrons may be limited by CR losses, where general formula 
(|29|) , that does not take into account the radiation reaction of accelerating particles, may not be 
applicable. 



After substituting expression ( |3Q ) for \p e \ and the above expression for $ into formula (|2 
and performing integration over £ and <fi, we get 

Vim(r?)| cosxral = - «) (l - ^) E rot . (34) 

Formally, in the above formula we should put r] — > oo to calculate the maximum power in the 
primary beam. Thus, for the maximum power of the primary beam we can write (see MH97, eq. 
[76]) 

3 

-f'prirrijmax = K )-^rot- (35) 

To estimate the pulsar bolometric photon luminosity it is reasonable to assume that 
L 7 0.5L pr i mjmax , where £ pr im,max is given by formula (|35|). Then, using expression (Q) 
for k we can write 

•-'45 (-. „,r'45 



L-y « 0.05-^1 (^1 - 0.15-^J ,B rot . (36) 

Thus, the energetics of CR photons generated by accelerating electrons above the pulsar PC is 
proportional to the pulsar spin-down luminosity, and, according to formula (|36|), the maximum 
efficiency of conversion of pulsar spin-down power into the high-energy quanta may amount to 
10%. Note that this estimate of efficiency implies that only half of the power of accelerating 
electrons gets consumed by 7-ray photons. 

In Figure 4 we plot the predicted high-energy luminosity as a function of spin-down luminosity. 
In the CR pair regime, where eq. (|3^) applies, we have plotted the luminosity for the screened 
unsaturated case (see top eq.p2]|). Below the CR death line, E rot $10 ergs -1 , the luminosity 
for the unscreened case of equation (|3~2"|) applies. Also plotted are the luminosities of the pulsars 
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with detected high-energy emission and their predicted luminosities. The detected high-energy 
pulsars are all above the CR death line, although Geminga and PSR B1055-52 are just above the 
line. The detected high energy pulsars seem to follow the predicted relationship, L 7 oc E rot . We 
predict that this relationship will break to L 7 oc E rot at E rot % 10 34 ergs -1 , so that older pulsars 
will have lower predicted luminosities than what would be predicted by an extrapolation of the 
trend seen in the younger pulsars. Non-thermal high-energy emission has in fact not been detected 
from older nearby pulsars such as PSR B1929+10 and PSR B0950+08, although a thermal 
emission component has been detected from PSR B1929+10 (Wang & Halpern 1997) which may 
be due to PC heating (see HM01). EGRET upper limits for pulsed emission from these pulsars 
(Thompson et al. 1994) are plotted in Figure 4 and lie above the predicted luminosities, but 
not by much. GLAST should be able to detect 7-ray emission from these older pulsars and test 
the predicted L 7 oc E ro t dependence and its location. Note that the break value of E TO t in the 
L y (E TOt ) dependence depicted in Figure 4, can be estimated by equating the screened expression 
for L 7 «0.5L p ^(£ rot ) given by formula ( |32| ) at L 7jmax m 0.5L pr i mjmax , where -L pr im,max is given 
by formula p5|), 



5 x W 1S E^P-V* = 0.04 E Iot , (37) 



which gives 



^rot,break = 1.4 x 10 34 P^ 7 B 12 2/7 erg • s" 1 . (38) 
For P 0.1 s, B\2 ~ 4, which is shown in Figure 4, -Erot.break ~ 6.6 x 10 33 erg-s -1 . Because of the 

—2/7 

B 12 dependence -E ro t,break should generally occur at higher E TOt values for ms pulsars. 



There are a number of ms pulsars which, according to equation (32), should have observable 
high-energy emission. Several of these, including PSR J0437-4715, PSR J0030+0451, PSR 
J1824-2452 and PSR J0218+4232, have non-thermal pulsed X-ray emission components but have 
not, with the possible exception of PSR J02 18+4232 (see Kuiper et al. 2000), been detected 
in the 7-ray band. It is important to remember, however, that the predicted luminosity L 7 is 
a bolometric luminosity. The luminosity in a particular band is very dependent on the actual 
emission spectrum. The vF v spectra of the bright 7-ray pulsars all peak in the 7-ray or hard 
X-ray bands with typical high-energy turnovers around several GeV. The high-energy spectrum of 
the ms pulsars may be quite different. The bulk of the ms pulsars are near or below the CR death 
line, so that their accelerating electric field is unscreened. The accelerating primary electrons will 
reach an energy where CR losses are compensated by the acceleration energy gain (see also Luo et 
al. 2000), 

7cRR =1.8xl0 7 S 1 1 2 /4 4/ 1 4 5 P- 1 / 4 . (39) 

The CR emission spectrum of these electrons will be quite hard (photon index —2/3) and will not 
be cutoff by magnetic pair production at an energy of a few GeV, but will extend to the natural 
cutoff of the CR spectrum at 

^ - 3 h 7 - -2.3xlO^^-V /4 , (40) 



2 mc 2 p c V P 
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where p c is the radius of curvature of a dipole field line. The CR vF u spectrum of the ms pulsars 
is therefore expected to peak at 50-100 GeV energies. The numerical model spectra for millisecond 
pulsars of Bulik et al. (2000) have also shown this result. The high energy curvature radiation 
from ms pulsars therefore falls in an energy band that has been above that of satellite detectors like 
EGRET and below that of air Cherenkov detectors. There could be a second spectral component 
of almost comparable total luminosity due to the synchrotron pair cascade from ICS pairs, but 
the vF v spectrum of this component would peak at much lower energies. There is also a question 
of whether two-photon pairs could screen the accelerating field and limit the voltage to very low 
values. Our preliminary investigations indicate that screening by two-photon pairs occurs only for 
NS surface temperatures high enough (Tq ~ 4) that two-photon pairs from returning positrons 
is likely to disrupt the acceleration (see discussion at end of Section 2.3.3). This issue will be a 
subject for future, more detailed investigation. 

4. Summary and conclusions 

In this paper we have outlined the status of two basic regimes of primary particle acceleration 
above a pulsar PC in the general-relativistic version of the space-charge limited flow model: the 
regime of acceleration with subsequent electron-positron pair formation by CR and screening 
of the electric field, and the regime of unscreened acceleration with ICS pairs and the emission 
of energetic CR photons. For the pulsar physical parameters we discussed in this paper such 
as the pulsar spin period, NS surface magnetic field strength, and PC temperature, the pair 
production may involve both the one-photon (magnetic) and two-photon (mostly 7-ray ICS 
photons on thermal X-ray photons) mechanisms, with high-energy photons generated via CR and 
ICS processes. In HM02 we began analyzing the onset of pair formation and addressed how the 
corresponding pair-formation criterion transforms into the pulsar theoretical death lines. HM02 
presented the death line calculations for the canonical NS model (mass of 1.4 Mq, radius of 10 
km, and moment of inertia I = 10 45 g • cm 2 ) and for one-photon pair formation mechanism only. 
Here we extended our study of pulsar theoretical death lines to include the effect of NS mass and 
radius and, for the ms pulsars with relatively hot PCs (T ~ 3 — 4 x 10 6 K), to incorporate the 
mechanism of two-photon pair formation. In the present study we demonstrate that the effect of 
NS mass is important for the death line calculations in ms pulsars. In fact, if the ms pulsars are 
more massive post-accreting NSs spun-up in low-mass binary systems, then the phase space they 
occupy in the P — P diagram would be consistent with our theoretical death lines calculated for 
NSs with masses more than 1.4 Mq. 

For the regime of acceleration without CR pairs we calculate the 7-ray luminosity of the 
primary beam as a function of pulsar spin-down power. It is important that the form of the 
theoretical L 7 — E rot dependence is determined by the regime of primary acceleration. For 
example, for the regime of acceleration accompanied by pair formation capable of screening the 
accelerating field L 7 oc E r ^ t , whereas for the unscreened regime of acceleration L 7 oc E TOt . We 
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discuss our theoretical L 1 — E rot plots in the light of currently available pulsar 7-ray data, and 

predict that for pulsars in the regime of unscreened acceleration (with relatively low values of 

■ 1/2 

-E-rot) L 1 should turn down from the E r ^ t dependence (see Figure 4). 

The main conclusions of our study can be summarized in the following way 

• We revised the prescription for the derivation of pulsar theoretical death lines to include 
the effect of variation (by more than an order of magnitude) of pair-formation altitude with 
pulsar parameters P and P. 

• Pulsar theoretical death lines are strongly affected by the EOS of a NS, with the onset of 
pair formation facilitated in more compact NSs. 

• Theoretical death lines for ms pulsars produced for massive NS models are in a good 
agreement with the empirical death line and with the hypothesis that ms pulsars are 
descendants of accreting NS in low-mass binaries. 

• The ms pulsars with relatively hot PCs (temperature ~ 3 — 4 x 10 6 K) may allow the 
occurrence of two-photon pair production, which significantly facilitates pair formation and 
may also move down the theoretical death lines in the ms range. 

• We derive a theoretical expression for L 7 and illustrate that it scales simply as E TOt or 

"1/2 

E T ^ t , respectively, for the unscreened regime of acceleration of the primary beam or for the 
acceleration implying pair formation and screening of the accelerating electric field. We 
predict that the break in L 7 — E xot dependence, attributed to the transition from the regime 
of acceleration with pair screening to the unscreened regime, might be seen in observational 
data. 

We acknowledge support from the NASA Astrophysics Theory Program. 

REFERENCES 

Arons, J. & Scharlemann, E. T. 1979, ApJ, 231, 854. 
Bulik, T., Dyks, J. & Rudak, B. 2000, MNRAS, 317, 97. 
Chen, K., & Ruderman, M. 1993, ApJ, 402, 264 (CR93). 
Cheng, K. S., Ho, C. & Ruderman, M. A. 1986, ApJ, 300,500. 
Daugherty, J. K. & Harding, A. K. 1996, ApJ, 458, 278. 
Friedman, J. L., Ipser, J. R., & Parker, L. 1986, ApJ, 304, 115. 

Glendenning, N. 1997, Nuclear and Particle Physics of Compact Stars (Berlin: Springer). 



-18- 



Harding, A. K., & Muslimov, A. G. 1998, ApJ, 508, 328 (HM98). 
Harding, A. K., & Muslimov, A. G. 2001, ApJ, 556, 987 (HM01). 
Harding, A. K., & Muslimov, A. G. 2002, ApJ, 568, 000 (HM02). 
Hibschman, J. A., & Arons, J. 2001, ApJ, 554, 624 (HA01) 
Kapoor, R. S., & Shukre, C. S.2001, A&A, 375, 405 

Kozlenkov, A. A., & Mitrofanov, I. G. 1986, Sov. Phys.-JETP Lett., 64, 1173. 
Kuiper, L., et al. 2000, A&A, 359, 615. 

Luo, Q., Shibata, S. & Melrose, D. B. 2000, MNRAS, 318, 943. 

Manchester, R. N. et al. 2001, MNRAS, 328, 17. 

Muslimov, A. G., & Tsygan, A. I. 1992, MNRAS, 255, 61 (MT92). 

Muslimov, A. G.> k Harding, A. K. 1997, ApJ, 485, 735 (MH97). 

Romani, R. W., 1996, ApJ, 470, 469. 

Ruderman, M., & Sutherland, P. G. 1975, ApJ, 196, 51. 

Sturrock, P. 1971, ApJ, 164, 529. 

Svensson, R. 1982, ApJ, 258, 335. 

Thompson, D. J., et al. 1994, ApJ, 436, 229. 

Thompson, D. J., 2001, in High-Energy Gamma-Ray Astronomy, ed. F. A. Aharonian & H. J. 
Volk, (AIP, New York), p. 103. 

Umeda, H., Shibazaki, N., Nomoto, K., & Tsuruta, S. 1993, ApJ, 408, 186. 

Usov, V. V. 2002, |astro-ph/0204402 . 

Wang, F. Y.-H., & Halpern, J. P. 1997, ApJ, 482, L159. 

Xu, R. X., Qiao, G. J., & Zhang, B. 1999, ApJ, 522, L109 

Zhang, L., & Cheng, K. S. 1997, ApJ, 487, 370. 

Zhang, B. 2001, ApJ, 562, L59. 

Zhang, B., Harding, A. K., & Muslimov, A. G. 2000, ApJ, 531, L135. 
Zhang, B., & Qiao, G. J. 1998, A&A, 338, 62. 



This preprint was prepared with the AAS IATgX macros v4.0. 



-19- 

CR, ICS Pair Front 

10" 12 
10" 14 
r 10" 16 

C/D 

• rT 10 

10" 20 

10" 22 
10" 24 

10" 3 10~ 2 10~ 1 10° 10 1 

P(s) 

• 

Fig. 1. — Pair death lines in the pulsar P-P diagram for curvature radiation (unlabeled 
curve) and for inverse-Compton radiation for different NS equations of state (see text). The 
parameter £ = 9/6 pc , where 6 pc = ({IR/c) 1 / 2 , indicates the magnetic colatitude of the primary 
electron acceleration. Also shown are radio pulsars in the ATNF Pulsar Catalog (available at 
http: / /www.atnf .csiro.au / pulsar^ ) . 
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Fig. 2. — Pair death lines for inverse-Compton radiation in the pulsar P-P diagram of short-period 
pulsars from the ATNF Pulsar Catalog, for rotating NS models having different masses (see text). 
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Fig. 3. — One-photon and two-photon pair death lines for inverse-Compton radiation in the pulsar 
P-P diagram for different PC surface temperatures, T§ = T/10 6 K, as labeled. Also shown are 
radio pulsars in the ATNF Pulsar Catalog. 
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Fig. 4. — Predicted and observed high energy luminosity vs. spin-down luminosity. The solid 
curve is the theoretical prediction from first expression in eq. ( |32| ) and eq. (|36|). The solid circles 
are the luminosities of the detected 7-ray pulsars, Thompson (2001), derived from detected fluxes 
above 1 eV assuming a 1 sr. solid angle. The upper limits are for > 100 MeV from Thompson et 
al. (1994). The open triangles are predicted luminosities for the detected pulsars. 



